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Abstract 

We report the first observation of Shubnikov-de Haas (SdH) oscillations and quantized Hall 
resistance in the multilayered massless Dirac fermion system a-(BEDT-TTF)2l3 with tilted cones. 
Holes were injected into the thin crystal fixed on a polyethylene naphthalate (PEN) substrate by 
contact electrification. The detection of SdH oscillations whose phase was modified by Berry's 
phase 7r strongly suggested that the carrier doping was successful in this system. We succeeded in 
detecting the quantum Hall effect (QHE) with the steps which is the essence of two dimensional 
Dirac fermion systems. The number of effectively doped layers was examined to be two in this 
device. We reveal that the correlation between effective layers plays an important role in QHE. 

PACS numbers: 73.43.Fj, 72.80.Le 



1 



Since the massless Dirac fermion system was realized in graphene [1, 21, the Dirac particles 
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discovered in other materials, such as graphite |3|, |4J, bismuth |5j, Fe-based superconductors 
6], topological insulators (7], and organic conductors {si^], have also fascinated physicists. 
Among them, the organic conductor a-(BEDT-TTF) 2 l3 [l6| (Fig. 1(a)) under high pressure 
offers an ideal testing ground for the Dirac particles. In this paper, we proffer the novel 
carrier doping technique for the massless Dirac fermion system a-(BEDT-TTF)2i3 under 
high pressure. Using this technique, we first observed Shubnikov-de Haas (SdH) oscillations 
and quantum Hall effect (QHE) in this system. 

«-(BEDT-TTF)2l3 is the first bulk (multilayered) two-dimensional (2D) zero-gap conduc- 
tor with massless Dirac particles. One of the characteristic features of multilayered systems 
is seen in the inter-layer transport under a transverse magnetic field. In the magnetic field, 
the energy of Landau levels in a zero-gap structure is expressed as 
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where vp is the Fermi velocity, n is the Landau index, and B is the magnetic field strength. 
One important difference between zero-gap conductors and conventional conductors is the 
appearance of a (n=0) Landau level at zero energy 17J. This special Landau level is called 
the zero mode. The characteristic features of zero-mode Landau carriers are clearly seen 
in inter-layer transport {ill, fisj] . In this field configuration, the effect of the magnetic field 
appears only through the change in the zero-mode density because the interaction between 



the electrical current and the magnetic field is weak. Thus, remarkable negative (inter 
magnetoresistance proportional to the inverse magnetic field was discovered 111 . 
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addition, recent progress in specific heat [l9| or NMR 20] studies also revealed characteristic 
feature of bulk Dirac system. 

Another significant feature different from graphene is that the Dirac cones in the present 
system are tilted, which has been indicated by the tight-binding band calculation 



the first principles band calculation 
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14| and 



15l | . Hence, the Fermi velocity, vp, is highly anisotropic. 



The relativistic Landau level structure (Eq. (pQ)) can be corrected, replacing v-p with the 
average vp. Thus, the tilt of the Dirac cones contributes to the small value of vp, it was 
estimated to be approximately 3.5 x 10 4 ms _1 [3]. This value is approximately 1/30 times 
that of graphene. 

Thus, this material with layered structure and tilted Dirac cones belongs to a broader 
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FIG. 1: (color online), (a) Crystal structure of a-(BEDT-TTF)2l3 viewed from the a-axis. (b) 
Schematic diagram of this system. The thickness of the crystal measured with a step profiler was 
approximately 100 nm. (c) Schematic energy diagram of the present device. The energy spectra 
(Dirac cones) for each layers are drawn. 

category of 2D massless Dirac fermion system. However, the relativistic Landau level (Eq. 
(CTJ) which is the most important characteristic feature in the Dirac fermion system has not 
directly been obtained yet. Both SdH oscillations and QHE have not been observed yet in 
this system until now, because the Fermi level is always located at the Dirac point. Moreover, 
the multilayered structure makes control of the Fermi level by the field-effect-transistor 
(FET) method much more difficult than the case of graphene. In order to clarify the primitive 
properties of our Dirac fermion systems, however, the detection of SdH oscillations and QHE 
are crucial. In this work, we made a breakthrough in the detection of SdH and QHE in the 
multilayered massless Dirac fermion system a-(BEDT-TTF) 2 l3 under high pressure. Our 
idea is the following. 

According to our investigation of the Hall effect in a-(BEDT-TTF)2i3, the carrier density 
at low temperature is only 10 8 cm _2 /sheet 



101 ] . Yet, the carriers are not localized but mobile 
with high mobilities. Thus, by fixing a crystal on a substrate weakly negatively charged by 
contact electrification, the effects of hole doping can be detected in the transport. Indeed, we 
succeeded in detecting the hole doping effects on the magnetoresistance and the Hall effect 
by fixing a crystal onto a polyethylene naphthalate (PEN) substrate (Fig. 1(b)). Figure 1(c) 
is the prospective energy diagram of this device. Holes should be injected into a few layers 
(pairs of BEDT-TTF molecular layers and 1% anion layers; vide infra). From this study, we 
anticipate that the surface charge produced by contact electrification of the PEN substrate 
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under high pressure is higher than 10 3 Cm 2 at low temperatures. 



Single crystals of a-(BEDT-TTF)2l3 were synthesized by the electrolysis method 16]. 
The crystal with a thickness approximately 100 nm was tightly fixed on the PEN substrate 
by van der Waals forces and/or electrostatic forces, as shown in Fig. 1(b). The resistivities 
and the Hall resistivities at p = 1.7 GPa were measured in magnetic fields of up to 7 T at 
temperatures below 2 K. Experiments were conducted as follows. A sample to which eight 
electrical leads were attached was put in a Teflon capsule filled with a pressure medium 
(Idemitsu DN-oil 7373), and then the capsule was set in a clamp- type pressure cell made of 
the BeCu hard alloy. 

Figure 2(a) shows the magnetic field dependence of the resistance at temperature below 
1.7 K. We notice an oscillation in the resistance at magnetic fields above 2 T. The amplitude 
of the oscillation increases with decreasing temperature. One vital clue is that the resistance 
oscillates as a function of B" 1 , as shown in Fig. 2(b). We succeeded in detecting the SdH 
oscillations with a frequency of approximately 8.3 T (= Bq) in this system for the first time. 

We demonstrate that the SdH oscillations in Figs. 2(a) and (b) originate from a 2D 
Dirac-type energy structure as follows. In the case of a Dirac fermion system, the circular 
orbit around the Dirac point in the magnetic field would yield Berry's phase 7r. This phase 
would modify the Landau quantization condition so that the zero-mode (n = 0) Landau level 
would always appear at the Dirac point. The effect of this phase is further probed in the 
semi-classical magneto-oscillation description, such as SdH oscillations. The SdH oscillation 
part of resistivity written by 

AR XX = A(B) cos [2ir(B /B + 1/2 + 7)] (2) 

acquires the phase factor 7=0 or 1/2 for normal electrons with Berry's phase and Dirac 
particles with Berry's phase 7r, respectively [2J. Here, A(B) is the SdH oscillation amplitude. 
We obtain phase factor 7 by plotting the values of B' 1 at the oscillation minima of AR XX 
(Fig. 2(b)) as a function of their number, as shown in Fig. 2(c). The linear extrapolated 
value close to 1/2 of the data to B~ l = determines the phase factor, 7=1/2. Thus, we 
identify the Dirac particles with a strong 2D nature from the angle dependence of SdH 
frequencies that obeys the 1/ cos(6 l ) law, as shown in Fig. 2(d). The polarity of the Hall 
resistance R xy , on the other hand, is positive in the magnetic field region in which SdH 
oscillations are observed, as shown in Figs. 3(a) and 4(a). Those findings strongly evidence 
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FIG. 2: (color online), (a) Magnetic field dependence of the resistance at temperatures below 
1.7 K. Magnetoresistance oscillations are observed. The presence of oscillations with temperature- 
independent frequencies proportional to B^ 1 shown in (b) strongly suggests that these are SdH 
oscillations, (b) SdH oscillations as a function of B^ 1 at 0.5 and 1.7 K. (c) The value of -B -1 for 
the SdH oscillation minima n m i n of AR XX (close circle) and for the filling factor v (open triangle), 
(d) Angle dependence of SdH frequencies. 

the following. This material fixed on the PEN substrate under high pressure is a 2D mass- 
less Dirac fermion system. Though the Dirac cone is tilted, the Landau level structure is 
expressed as Eq. Moreover, the hole doping in this system is a success. 

Here, we discuss the hole doping level of this device. In the device that uses a multilayered 
material, carriers are expected to inject into a few layers. According to Horowitz, et al, the 
doped carrier density n(z) in the multilayered materials is written by 



n(z) oc (N Z + LY 



(3) 



from the Poisson's equation, where N z is a number of layers from an interface and L can 
be regarded as the effective layers (thickness) 2l|. The simple SdH signals shown in Fig. 
2(a) suggest that "the principal effective layer" is a single. If several pairs were operated, we 
would be able to obtain much more complex SdH signals because each doping level would 
be different (Eq. Performing the second-order differential of R xx , however, we find the 

SdH oscillations with two frequencies as shown in Fig. 3(a). One is the first observed SdH 
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FIG. 3: (color online), (a) Magnetic field dependence of R xx , R xy and — (d 2 R xx /dB 2 ) at 0.5 K. The 
arrows shows second Sdh oscillation minima of — (d 2 R xx /dB 2 ). (b) The value of B~ l for the SdH 
oscillation minima n m j n of —(d 2 R xx /dB 2 ). (c) The distribution of doped carrier density to layered 
direction. The line is the curve of Eq. ([3]) with L = 2. 

oscillation with Bq ~ 8.3 T and another one (which is denoted by arrows in Fig. 3(a) ) has 
a frequency of B ~ 20.4 T. The values of B~ x for SdH oscillation minima of —(d 2 R xx /dB 2 ) 
are shown in Fig. 3(b). The linear extrapolated value of another one close to 1/2 of data to 
B^ 1 = strongly indicate that the origin of this SdH oscillation is also Dirac particles. Thus, 
the holes are injected to more than two layers, and two layers from an interface emitted the 
SdH oscillations with different frequencies respectively. 

Let us examine the distribution of doped-holes density to layered direction. The 
carrier density of each layers are estimated to be approximately 19.7 x 10 n cm -2 and 8.3 x 
10 11 cm -2 from the relationship = Bog/<f)Q, where g=4 is the four-fold degeneracy and 
0o=4. 14x 10~ 15 Tin 2 is flux quantum. The Fermi surface areas are approximately 0.05% and 
0.02% of the first Brillouin zone under ambient pressure. Judging the distribution of the 
carrier density from Eq. ([3]), of the first layer from an interface is higher one (19.7 x 10 11 
cm -2 ). Thus, Fig. 3(c) is the distribution of doped carrier density to layered direction in 
this device. The energy diagram such as Fig. 1(c) is expected. When we assume that the 
Eq. (j3J) can be adapted for our system, the best fit of this equation with L = 2 indicates 
that the number of effective layer is two. Note that primary (simple) SdH signals shown in 
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Fig. 2(a) are due to the Dirac particles of the second layer because the carrier mobility in 
the first layer is low. The lower limits of the carrier mobilities in the first and the second 
layers are roughly examined to be approximately 3.8 x 10 3 cm 2 V _1 s _1 and 10 4 cm 2 V _1 s _1 
from the early SdH signals (u c r = {[iB) 2 > 1), respectively. Non-smooth surface of PEN 
will give rise to the low mobility in the first layer. 

The detection of the SdH oscillations motivated us to investigate QHE in this system. 
We show the magnetic field dependence of the Hall resistance R xy at 0.5 K in Fig. 2(a). 
The two obvious plateaux observed at magnetic field around 3.5 and 5.5 T show that R xx 
minima are the hallmarks of QHE. We also find the anomalies at fields around 2 and 2.5 T. 

Note that R xx is not zero but shows minima at the plateaux of R xy , as shown in Fig. 
2. In general complete 2D system cases, R xx should be zero at the quantum Hall plateaux 
because a xx = 0. The present system, however, has a multilayered structure even though 
few layers, primary the second layer is operated. In this case, the other layers contribute to 
the conduction so that R xx shows finite values at the quantum Hall plateaux which blurs 
the SdH oscillations. Moreover, undoped layers also yields the temperature dependence of 
R xy shown in Fig. 4(a). In this situation, it is difficult to directly demonstrate that the 
plateaux are complete quantum Hall states. 

The following findings, however, indicate that this is an intrinsic QHE: (1) The carrier 
mobility of this material is expected to be higher than 10 5 cm 2 V~ 1 s~ 1 at low temperature. 
The lower limit of the carrier mobility for the operated layer in this device, on the other hand, 
is expected to be ~ 10 4 cm 2 V _1 s _1 because SdH oscillations are clearly detected around 1 
T from the second order differential of R xx as shown in Fig. 3(a). In graphene or oxide 
leterostructures, QHE has been observed in samples with mobilities lower than those values 
l|, 0|. (2) As mentioned below, the steps of the plateaux are the same as QHE in 2D Dirac 
fermion systems l|, |2j . Especially, the clear plateaux of the Hall conductivity shown in the 
inset of Fig. 4(a) strongly indicates that those are the hallmark of intrinsic QHE. (3) The 
breakdown of the plateaux due to the high current intensity is the characteristic feature of 



QHE 23M25|| . We verified that the plateaux were blurred with increasing current intensity, 



as shown in Fig. 4(b). Note that the data of R xy in the normal state are independent of the 
current intensity. In the quantum Hall state, Joule heat is locally generated in the vicinity 
of current electrode at the sample edge. Nonequilibrium electron distribution generated at 
the sample edge gives rise to energy dissipation so that the breakdown of QHE occurs. 
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FIG. 4: (color online), (a) Magnetic field dependence of the Hall resistance at temperature below 
1.7 K. Arrows indicate the quantum Hall plateaux or anomaly. The inset is v dependence of the 
Hall conductivity a xy = p xy /(p xx + Pxy)- The plateaux are clearly observed at \u\ = 6 and 10. 
(b) Hall resistance for several current intensities at 0.5 K. The insets show the regions from 2 to 
3 T and from 3 to 4 T. The high intensities of the currents blur the plateaux, (c) Magnetic field 
dependence of the Hall resistance divided by 6 and the value at 6 T. Solid lines show the filling 
factors for v =-6, -10, -14, and -18, respectively. 

Next, we examine the plateaux of R xy . In Dirac fermion systems, such as graphene, R xy 
quantization is in accordance with R~y = v ■ e 2 /h, where v = ±g(n + 1/2) with g = 4 is the 
quantized filling factor. A conspicuous effect on the Dirac fermion system is that the factor 
of half-integer exists. Thus, probes of the quantum Hall plateaux for \u\ — 2, 6, 10, 14, 18, 
• • • are expected. 

The multilayered structure with a few hole-doped layers, however, gives rise to the lack 
of validity of the examination of v from the values of R xy plateau. The values of the R xy 
plateau depend on the thickness. We cannot ignore the sum of Hall conductivities for the 
undoped layers. We can see this effect on R xy in magnetic fields below 1.5 T. In this field 
region, the polarity of R xy is negative and yet holes are injected. This is frequently observed 
behavior of the present material induced by doped electrons with ppm order that originate 
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from unstable I J anions [6j. 

Thus, we easily examine the filling factors as follows. First, we assume that this material 
is truly a 2D massless Dirac fermion system, the same as graphene. Second assumption is 
that the plateaux or anomaly of R xy shown in Figs. 3(a) or 4(a) are the QHE operated on 
the second layer. The contribution of the first layer is not small. In order to easily examine 
the filling factor, however, we ignore the Hall conductivity of the first layer here. Based on 
these assumptions, the filing factors v = -6, -10, -14, -18 are expected at minima of R xx from 
the SdH oscillations as shown in Fig. 2(c). For example, we can regard the filling factor of 
the plateau around 5.5 T as z/=-6. When the filling factor of the plateau around 5.5 T is 
fixed at z/=-6, we find those of the other plateaux or anomaly to be close to v =-10, -14, and 
-18, respectively, as shown in Fig. 4(c). These steps are the essence of QHE in 2D Dirac 
fermion systems. Furthermore, the plateaux of a xy at v — 6 and 10 shown in the inset of 
Fig. 4(a) encourage us that those are the essence. 

Lastly, we briefly mention the bilayer system in our device. In the above discussion, the 
Hall effect of the first layer was ignored. Holes are injected into the first layer, however, 
though the mobility is lower than that of the second layer. In the case that the first and 
the second layers do not depend on each other, the quantum Hall state with a relationship 
ndx/nci2 = v\jvii should be stable at the filling factor of v — V\ + U2, where Udi and n$i are 
the carrier densities and v\ and are the filling factors at the first and the second layers, 
respectively. For example, we obtain Udi/n^ ~ ^1/^2 with v\ = —14 and ^2 = ~ 6 at the 
plateau for B ~ 5.5 T. The filling factor of the data at 0.5 K shown in Figs. 3(a) or 4(a), 
however, lacks validity of this relationship v — v% + z/2, it is v — —11 (7^ —14 — 6). The 
consideration of undoped layers strengthens this inconsistency. This discrepancy suggests 
that the correlation between the first and the second layers plays an important role in 
the QHE. Recently, Morinari theoretically demonstrated that the correlation between layers 
narrows and blurs the plateau 26]. In the bilayer device with high carrier mobilities of GaAs, 
on the other hand, peculiar quantum Hall states associated with the coherence between layers 
was discovered 27j. Thus, the present device with higher carrier mobilities would provide a 
new type of quantum Hall state under higher magnetic fields. 

In conclusion, we provided crucial evidence that a-(BEDT-TTF) 2 l3 under high hydro- 
static pressure is an intrinsic zero-gap conductor with Dirac cone type energy structure. This 
is the first work that demonstrated the carrier injection effects on the transport phenomena 
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in the bulk zero-gap conductors. We achieved the hole doping into a-(BEDT-TTF) 2 l3 by 
fixing the thin crystal to the PEN substrate. The number of effectively doped layers was 
examined to be two. We were able to detect SdH oscillations with two frequency whose 
phase are modified by Berry's phase tt. The quantum Hall plateaux for z/=-6, -10, (-14, and 
-18) due to the primary effective layer was observed. These steps are the essence of QHE in 
2D Dirac fermion systems. Moreover, we revealed that the correlation between bilayer plays 
an important role in QHE. This device offers an ideal testing ground for the bilayer QHE of 
the massless Dirac fermion systems. 

This work was supported by Grants-in-Aid for Scientific Research (No. 22540379, 
No.22224006, and No. 23740270) from the Ministry of Education, Culture, Sports, Sci- 
ence and Technology, Japan. 
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